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The first-principles method and the evolutionary algorithm are used to identify stable high pressure phases 
of potassium azide (KN 3 ). It has been verified that the stable phase with space group 14/ mem below 22 GPa, 
which is consistent with the experimental result, will transform into the Cllm phase with pressure 
increasing. These two phases are insulator with anions. A metallic phase with P6/mmm symmetry is 
preferred above 40 GPa, and the N atoms in this structure form six-membered rings which are important 
for understanding the pressure effect on anions and phase transitions of KN 3 . Above the studied 
pressure (100 GPa), a polymerization of N 6 rings may be obtained as the result of the increasing 
compactness. 

Metal azides have drawn considerable attention for their interesting chemical and physical properties. 
Under external influences (impact, heat, irradiation, etc), they become unstable and decompose into 
metal and nitrogen 1 . Their practical applications include explosives, pure nitrogen sources and pho- 
tographic materials. Additionally, metal azides are structurally simple among solids that deflagrate or detonate, so 
they are potentially model systems for theories of the fast reactions 2 . Being different from the extensively studied 
alkali halides, metal azides are ionic compounds containing internal molecular structure which makes them the 
candidates for understanding the complex nature of chemical bonding. On the other hand, the polymerization of 
nitrogen may form a high- energy- density material because the transformation from the N-N single bond 
(160 kj/mol) to the N=N triple bond (954 kj/mol) is accompanied by a large energy release. The non-molecular 
nitrogen under high pressure was first predicted theoretically by A. K. McMahan et aV. A single-bonded cubic 
gauche form of N 2 (cg-N) was also successfully synthesized by M. I. Eremets et al. 4,5 . It is suggested that cg-N can 
be stabilized in compounds with other elements or by introducing impurities 4 . Recently, metal azides have been 
proposed to be a precursor in the formation of polymeric nitrogen. It could be expected that the anion will 
create polymeric single- covalent-bond networks more easily than diatomic nitrogen because the N3" anion is 
more weakly bonded than the diatomic triple-bonded nitrogen 6 . Alkali azides are one class of compound among 
metal azides. Pressure-induced phase transitions in alkali azides have been reported by both experimental 6 " 11 and 
theoretical investigations 12 " 16 . Previous study reported that sodium azide undergoes a set of phase transitions 6 ' 8 , 
and the N3" anions in NaN 3 transform to polymeric nitrogen net above 120 GPa 6 . Lithium azide, which is 
isostructural to the low-temperature phase of NaN 3 at ambient condition, is stable up to the pressure of 
60 GPa at room temperature 7 . Recent high-pressure study of CsN 3 up to 55.4 GPa reveals three phase transitions 
approximately at 0.5, 4.4, and 15.4 GPa 11 . Strikingly, theoretical studies have predicted polymerization of nitro- 
gen in LiN 3 14,15 and NaN 3 16 . Thus, a study of the high-pressure behavior of KN 3 would provide more insights into 
the mechanism of pressure-induced rearrangement of azide anions. It is helpful to investigate theoretically the 
pressure effect on potassium azide which might result in the formation of polymeric nitrogen. 

At ambient condition, potassium azide crystallizes in a body-centered tetragonal lattice with I4/mcm sym- 
metry, and the N3" ion in the lattice is a straight chain of three nitrogen atoms. Raman scattering up to 4.0 GPa 17 
and single-crystal x-ray diffraction up to 7.0 GPa 18 show that no phase transition is found at these measured 
pressures. Recently, a structural phase transition at 15.5 GPa is revealed by X-ray diffraction study 9,10 . A Raman 
scattering study of KN 3 up to 55.0 GPa suggests that a first-order phase transition starts at 13.6 GPa and 
completes at 32.2 GPa 10 . However, a detailed structure of solid KN 3 under high pressure has been unsolved 
yet. In this paper, we employ the first-principles study to understand the phase transitions and properties of KN 3 
under high pressure. 
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In the present work, KN 3 is compressed up to 100 GPa. The phase 
transition from 14/mcm to C2/m is determined at a pressure about 
22 GPa. In addition, we have defined another high-pressure phase 
with P6/mmm symmetry in which the N3" ions transform to N 6 
rings, and the hexagonal phase is metallic. We have just noticed 
the discussion on the P6/mmm phase by Li et at. recently 19 . The three 
phases of this crystal are stable against decomposition of KN 3 into 
3 

K+ -N 2 under pressure up to 100 GPa. The results provide an 
insight into the formation of polymeric nitrogen in metal azides. 

Computational details 

To find the lowest energy structures of KN 3 under extreme condi- 
tions, the USPEX code based on the evolutionary algorithm 20,21 is 
employed for the search where the VASP code 22 is used as an external 
ab initio code for the underlying structural optimizations. 
Evolutionary variable- cell structure prediction simulations are per- 
formed at 20, 60, and 100 GPa with 1, 2, and 4 KN 3 formula units per 
unit cell. The first generation is produced randomly. All newly gen- 
erated structures are relaxed at constant pressure and ranked by their 
enthalpy value. The lowest-enthalpy 60% structures of each genera- 
tion are used to produce the next generation through heredity (70%), 
atomic permutation (10%), atom position mutation (10%) and lattice 
mutation (10%). The structures are relaxed by using density- func- 
tional theory with the Perdew-Burke-Ernzerhof (PBE) exchange cor- 
relation functional 23 , and projector augmented wave (PAW) 
method 24 is also adopted. The 3s 2 3p 6 4s 1 for K and 2s 2 2p 3 for N are 
treated as valence electrons. Having selected the lower enthalpy 
structures, we recalculate their enthalpies with increased accuracy 
between 0 and 100 GPa. An energy cutoff of 520 eV is used for the 
plane-wave basis sets, and the total energy is converged to 1.0 X 
10" 6 eV in the self-consistent loop. The atomic positions, lattice 
parameters, and cell volume are fully optimized by using a conjug- 
ate-gradient algorithm. The iterative relaxation of the atomic posi- 
tions stop when all forces are smaller than 0.001 eV/A, and the total 
stress tensor is reduced to the order of 0.01 GPa. The phonon calcu- 
lations are performed using the Quantum ESPRESSO code based on 
density functional perturbation theory 25 , where vanderbilt ultrasoft 
pseudopotentials are used for K and N. Before the application, the 
pseudopotentials are tested by comparing the relaxed structural 
parameters and electronic structure with the results obtained from 
VASP code. All cases concerned are zero -temperature ground state. 

Results and discussion 

The analysis of the predicted structures gives us a list of candidate 
structure with space groups 14/mcm, C2/m, P2 l5 Cmc2 l5 P6/mmm, 
Pl y and Imma, which are depicted in Fig. 1. In order to exam the 
difference between theoretical results and experimental data, the 
theoretical lattice constants and unit cell volume are compared with 
experimental I4i 'mem structure at ambient condition (Table I). The 
LDA calculations underestimate a, c, and V by 3.3%, 5.2%, and 
11.4%, whereas GGA calculations overestimate a, c, and V by 1.5%, 
2.0%, and 5.2%, respectively. The GGA results are closer to the 
experimental results than the LDA ones. Therefore the GGA 
exchange-correlation functional is adopted for the further calcula- 
tions. The C2/m structure detected from the simulations has the 
same space group as in LiN 3 7 , NaN 3 26 and CsN 3 n . The 14/mcm, 
C2/m, P2 X and Cmc2 l structures possess linear molecular N3" anions, 
while the structures with Pi and Imma symmetries display the char- 
acteristic of the N-atom chains and are radically different from the 
structures containing N3" anions (see Fig. 1 and supplementary 
Table SI). In variable-cell simulations at 100 GPa, the P6/mmm 
structure is clearly identified as the most stable one and it is assigned 
to a new structure type for alkali metal azides. Interestingly, this is a 
layer-like structure with six N atoms forming a planar N 6 ring. 
Additionally, polynitrogen molecules as clean high energy density 
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Figure 1 | The structures computed for KN 3 in the pressure range 
0-100 GPa. The purple spheres are potassium atoms and grey are nitrogen 
atoms. The 2X1X2 supercell of the PI structure and the 2X1X1 
supercell of the Imma structure are viewed along b-axis. 

materials have been extensively explored for several years 27 " 29 . A lot 
of theoretical work on N 6 has shown that the planar hexagonal (D 6h ) 
ring is not minimum 30 " 32 . However, the six-membered nitrogen rings 
can be stabilized by coordinate covalent to oxygen 33 and by the 
incorporation of metal atom in metal-N 6 molecules 34 . We now 
obtain the N 6 rings in bulk materials. It reveals that extra atoms play 
an important role in stabilizing N 6 ring. More recently, a P6/m struc- 
ture containing N 6 ring has been also predicted in LiN 3 13 " 15 and 
NaN 3 16 . Thus, the P6/m structure is considered in KN 3 . The results 
show that the P6/m structure becomes P6/mmm phase after structure 
optimization. Parameters describing the C2/m and P6/mmm struc- 
tures are listed in Table II. 

The enthalpies of the most energetically competitive structures are 
compared over the pressure range 0-100 GPa as shown in Fig. 2. The 
most stable structure is a tetragonal phase with 14/mcm symmetry 
from ambient pressure up to 22 GPa, which is then replaced by a 
lower- enthalpy C2/m structure. Actually, at 13.6 GPa, Raman spec- 
tra 10 have identified a phase transition that completes at 32.2 GPa, 
which is in agreement with the C2/m phase being stable above 
22 GPa. The I4/mcm-C2/m sequence has also been observed in 
CsN 3 u . During I4/mcm-C2/m transition, the tetragonal lattice is 
distorted under pressure, and the orientation of N3" anions changes 
between the two K layers. The N3" anion in the C2/m structure is 
parallel to one another. Moreover, compression induces a symmetry 
reduction and N3" anions are still in a molecular state in this trans- 
ition, which are consistent with the experimental results 10 . Above 
40 GPa, a hexagonal structure with P6/mmm symmetry is favored 
over other structures and remains the lowest-enthalpy phase up to 
100 GPa. We find that the crystal structures of KN 3 containing the 
N3" ions are energetically favorable at lower pressure, while at higher 
pressure, there is a tendency to the N chains or rings. Although 
pressure induces the rearrangement of azide ions, the formation of 
N-N single bond needs further compression. The dependence of 
volume on pressure is shown in the inset of Fig. 2. The volume 
reductions of —2.3% and —9.5% are found for the I4/mcm-C2/m 



Table 1 Theoretical lattice constants and unit cell volumes V at 
ambient pressure compared with experimental data for 14/mcm 
structure 






Present work 


Parameter 


Experiment (ref. 9) 


LDA GGA 


a = b(A) 

c(A) 

V(A 3 ) 


6.1 1094 
7.09755 
265.043 


5.90883 6.20556 
6.72696 7.24091 
234.867 278.841 
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Table II Optimized structural parameters of C2/ 
tions 


77 phase at 30 GPa and P6/mmm phase 


at 1 00 GPa from the first-principles calcula- 


Pressure(GPa) 


Space group 


Lattice parameters! A, °) 


Atomic coordinates(fractional) 


30 
100 


C2/m 
P6/ mmm 


a = 4.293, k = 4.360, c = 4.767 
oc = 90 t p= 107.969,7 = 90 

a= b = 5.376, c= 2.366 
a = p = 90, y = 120 


K 2a (0.0000, 0.0000, 0.0000) 
N 2d (0.0000, 0.5000, 0.5000) 
N 4/(0.6272, 0.0000, 0.3151) 
K 2d (0.6667, 0.3333, 0.5000) 
N 6/ (0.0000, 0.2396, 0.0000) 



and C2/m-P6/mmm transitions, respectively. The discontinuous 
change in volume indicates that the two phase transitions are first 
order. As the transformation from C2/m to P6/mmm involves the 
forming of N-N bonds and is reconstructive, there is a large kinetic 
barrier. In addition, the possibility of decomposition formula of 
3 

KN 3 = K+-N 2 is checked by the enthalpies of decomposition, 

where we consider Fm3m, Pnma, MJamd, and Cmca structures 
for K 35 ' 36 and the cubic gauche {I2 Y 3) structure for the N 2 4 ' 37 . It turns 

3 

out that the enthalpies of K + - N 2 are much higher than that of KN 3 

at the concerned pressure range (>3.2 eV). Thus the KN 3 crystal 
keeps stable against decomposition over the 0-100 GPa pressure 
range. 

Since the high-pressure P6/mmm phase is novel for alkali metal 
azides, it is essential to investigate the properties of P6/mmm struc- 
ture in detail. The ratios of a/a 0 and c/c 0 for the P6/mmm structure 
seen in Fig. 3a exhibit that the compression is anisotropic with the 
reduction of lattice parameter a by 13.97% over the pressure studied, 
while lattice parameter c decreases by 26.73%. This implies that it is 
more compressive along interlayer direction than intralayer dir- 
ection for the layer-like P6/mmm phase. The N 6 ring with D 6h sym- 
metry possesses benzene-like characteristics. The charge density of 
N 6 ring in Fig. 3c suggests that each nitrogen atom forms two a bonds 
with its two nearest neighbors by overlapping two sp 2 hybrid orbitals, 



which is similar to the C-C bonds of C 6 H 6 . Additionally, for covalent 
bonds, there is a general trend in that the shorter the bond length, the 
stronger the bond strength. As shown in Fig. 3b, the N-N bond length 
of N 6 ring is longer than that of N=N double bond (1.25 A) but 
shorter than that of N-N single bond (1.45 A). This means that the 
a bond of N 6 ring is weaker than the N=N double bond but stronger 
than the N-N single bond. In general, pressure induces the destabi- 
lization of intramolecular bonds. The structure with N 2 and N 6 units 
has been reported around 60 GPa by first-principles simulations for 
pure nitrogen 38 . Recently, a diamondoid structure of the polymeric 
nitrogen is predicted above 263 GPa 39 . These findings also provide a 
basis for understanding the high-pressure behavior of nitrogen- 
related materials. As pressure increases, there is a competition 
among the a bond, ionic bond, and van der waals. Thus, it is probably 
fair to say that N 6 rings in P6/mmm phase will transform to poly- 
meric nitrogen networks under further compression. 

The ambient-pressure phase 14/mcm is insulating with a calcu- 
lated band gap of 4.2 eV at 0 GPa (see supplementary Fig. SI), and 
the C2/m structure is also found to be an insulator with a band gap of 
4.1 eV at 30 GPa (see supplementary Fig. S2). On further compres- 
sion, KN 3 transforms to a metallic state with P6lmmm symmetry at 
40 GPa. Besides, the PI and Imma structures are metastable above 
60 GPa. The Imma structure is metal (see supplementary Fig. S3), 
and the Pi structure is a semiconductor with narrow band gap 
(0.78 eV, see supplementary Fig. S4). Fig. 4 presents the band struc- 
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Figure 2 | Calculated enthalpies per KN 3 unit as the function of pressure. The enthalpies are referenced to that of Cllm. Inset: the pressure dependence of 
volume for Wmcm, Cllm and P6/mmm phases of KN 3 . 
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Figure 3 | (a) Variation of normalized lattice parameters of P6/mmm structure with pressure. a 0 and c 0 are the lattice parameters of the equilibrium 
volume structure at zero pressure, (b) Pressure dependence of the N-N bond length of N 6 ring, (c) Total charge density plotted in the (001) plane for 
P6/mmm phase at 100 GPa. 



ture and the density of states of P6/mmm phase at 100 GPa. It reveals 
that P6/mmm structure is a weak metal with a small density of states 
at the Fermi level (0.07 states/eV/cell) which are associated with the 
N-2p electrons. The bands across the Fermi level are highly dispers- 
ive along the c-axis (G — » A, H —> K, and M —> L). In addition, the 
Bader method 40 is chosen to analyze the charge transfer as imple- 
mented in the algorithm developed by Henkelman et al. 41 . The Bader 
charges of the P6/mmm phase are around +0.73 and —0.24 for 
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Figure 4 | Electronic band structure and projected density of states 
(PDOS) for P6/mmm phase at 100 GPa. 



potassium and nitrogen, respectively, which suggests that the high- 
pressure P6/mmm phase has ionic characteristics for K-N chemical 
bonds. It means that the N 6 anion in P6/mmm phase has nearly 8 n- 
electrons. Six P z orbitals form three bonding n orbitals and three 
antibonding tt* orbitals. Thus, the n* orbitals are partially occupied 
by two electrons, accompanied by two conduction bands crossing the 
Fermi level. Furthermore, the dynamical stability of the P6/mmm 
structure is established from the phonon calculations. As shown in 
Fig. 5, the absence of any imaginary frequency confirms the stability 
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Figure 5 | Phonon dispersion curve and phonon density of states 
(PHDOS) for P6lmmm phase at 100 GPa. 
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of the P6/mmm phase. The lower bands, ranging up to 620 cm" 1 , are 
formed by a significant mixing of K and N vibrations, while the 
higher bands are mostly attributed due to N-atom. 

Conclusion 

In summary, an evolutionary algorithm in conjunction with first- 
principles electronic structure computations has been used to predict 
the stable high-pressure phases of potassium azide. Our calculations 
indicate that the experimental 14/ mem phase of KN 3 transforms to 
Cllm structure at 22 GPa and then to a hexagonal P6/mmm struc- 
ture at 40 GPa. The planar N 6 ring is formed in the metallic P6/mmm 
phase. This phase of KN 3 is dynamically stable. Further compression 
could lead to the polymerization of N 6 rings, which can be used as a 
potential high-energy-density material. 
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